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1. INTRODUCTION 

The problem of the ionization in the spiral galaxies, including our own, the Milky Way is 
not solved at the moment (for survey see e.g. [1]) although there are some indications that 
could lead towards solution [2, 3, 4], 

We assume the Sofue’s model of the rotation curves of spiral galaxies [5] based upon 
the observation of CO line emission and combined with HI and optical observations. This 
model has four mass components: nuclear mass component, a central bulge, a disk, and a 
massive halo. One can add the hfth component for the Milky Way [5] at the nucleus e.g. 
[6,7]. As a characteristic of the ionization we study the number density of free electrons 
in the interstellar medium throughout the Galaxy, Ug, using the Taylor and Cordes model 
(TC93) [8] as well as its updated version (LC98) [9]. 

2. THE CASE OF THE MILKY WAY 

We examine our Galaxy - Milky Way, in many respects a typical spiral galaxy, in more 
detail. In the textbooks one usually encounters the fact that in the disk of a spiral there 
are H II regions that are ionized by the radiation from young and massive stars, associa¬ 
tions of OB stars [10]. Although the average properties of the ionized gas are rather well 
known the source of this ionization is still the subject of the discussion of the astronomers, 
astrophysicists, and even cosmologists. Following Sciama [1] one can state the following 
important problems. 

• The scale height (dehned as the column density on one side of the galactic plane 
divided by the volume density in the plane) of free electron component of the inter¬ 



stellar matter (ISM) is hg 


where N is the the total column 


density {N ~ 6.8 x 10^^ cm“^ [11]) and (ne)o = 0.033 cm“^ is the mean electron 
density at the midplane. However, the sources such as O stars or supernovae have a 
much smaller scale height (hg ~ 100 pc). 

• As pointed out by Reynolds [12] the nature and origin of the ionized gas about the 
galactic disk is not well understood. The power that is required to maintain the 
diffuse, ionized gas ranges from minimum of 5.3 x 10“® ergs s“^ cm“^ (if all the gas 
has T ~ 10'^ K) to 6.6 x 10“^ ergs s“^ cm“^ (T ~ 10^ K). As stressed by [1] among 
all known sources of power only Lyc (Lyman continuum) photons from O stars and 
energy from supernovae have or exceed the required power. 
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• Therefore, since the interstellar medinm is highly opaqne to hydrogen-ionizing radia¬ 
tion, the following qnestion remains: how can this radiation travel hnndreds of parsecs 
from the parent O stars in order to prodnce the diffnse ionized gas? (e.g. [13]). Care- 
fnl stndy of the interstellar hydrogen towards different pnlsars snggests that colnmn 
densities of free electrons along two line segments cannot be acconnted for by H II 
regions aronnd B stars or hot white dwarf stars. Reynolds [12] gives three possibil¬ 
ities for the explanation of the existence of the ionized gas: (i) morphology for the 
interstellar H I is very different from that nsnally depicted, (ii) Lyc Inminosities for 
early B or hot white dwarf stars are more than an order of magnitnde larger that the 
cnrrently accepted valnes and (hi) there exists an additional, yet nnrecognized sonrce 
of ionization within the Galactic disk. 

• The mean electron density in opaqne interclond regions within a few hnndred parsecs 
of the Snn as f un ction of a Galactocentric radins is approximately constant [8]. 

TG93 modeled the electron density of a certain Galactic location as the following snm 

(the hfth component was added in LG98): 

ne{x,y,z) ^ nigi{r)sech‘^{z/hi) + n2g2{r)sedi^{z/h2) + 

UaSech^ {z / ha) E fjgair, Sj) -F nagaiu) + nGc 9 Gc{r)hGc{z) (1) 

The detailed description of eacti component is given in TG93 and LG98 and becanse of 
space limitation we jnst give few remarks: r is the Galactocentric distance projected onto 
the plane and is eqnal r = {x^ -|- the snm goes over fonr spiral arms, f = 1...5 

denotes the density in different regions, /^ , f = 1...4 are scale factors, gi ^ i — 1.. .b are 
fnnctions of position, hi, i = 1.. .4 are scale heights and 2 : is the height above the galactic 
plane. Electron densities as a fnnction of R are presented in TG93 (in their Fig. 3). 

The hfth component becomes dominant np to a Galactocentric distance of 0.6 kpc. 
In this area, nnmerons point sonrces of ionization are responsible for the total Galactic 
center (GG) component contribntion, especially within the nnclear bnlge (r < 250 pc). 

At a dynamical center of the Galaxy, SgrA* - a nonthermal synchrotron radio-sonrce 
coinciding with the estimated 2.6 x 10 ®Mq black hole [6], lies snrronnded with thermal 
plasma - SgrA West, ionized by a clnster of few tens of yonng, hot (Tgff = 3.5 x 10^ 
K), Inminons O stars within a central (r = 1 pc) cavity [14]. SgrA West consists of a 
minispiral and an extended component with n|P ~ 10^ cm“^ and ~ 10^ cm“^ [15]. 
Some local featnres of ionized gas snch as the Bnllet (4” northwest of SgrA*) [16] and the 
Sickle [17] contribnte to the total electron density near GG with 3 x 10^ cm“^ and 150 
cm“3 respectively. 

At a distance of 30 pc behind the GG, lays SgrA East - the probable snpernova 
remnant colliding with a molecnlar clond at 200 pc from GG, contribnting with Ue = 
6 cm“^ electron density. It conld be heated by O stars formed as a resnlt of collision [19]. 

Another large ionized area near GG is a Sgr B2 - the largest Inminons H Il/molecnlar 
clond region at 100 pc from GG snrronnding the circnmnnclear disc (2-12 pc from GG) in 
which the SgrA West is embedded. SgrB2 conld be interpreted as an extended accretion 
disk aronnd a central dormant black hole [20]. 

If we now wish to stndy the ionization in onter parts of the Galaxy we will enconnter 
the lack of star forming regions in the ontskirts of spiral galaxies [21]. From the radio ob¬ 
servations (21 cm) it was noticed that HI disks in at least two galaxies (M33 and NGG3198) 
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have sharp edges [1]. HI column density drops from few times 10^^ cm“^ to few times 10^® 
cm“^ within 1 to 2 kpc. When the incident photon flux F is monochromatic and the ion¬ 
ized gas is smoothly distributed (with density n and thickness b) the edge will occur when 
jy ^ F. where a is the recombination coefficient. If one wishes to estimate the intergalac- 
tic hydrogen-ionizing flux F, one can measure Ha surface brightness of clouds opaque to F 
exposed to this flux - each hydrogen recombination produces 0.46 Ha photons on average 
(e.g. [1]). After extensive studies an upper limit for ionizing flux is currently thought to be: 
Fext < 6 X 10^ cm“^ s“^ [1]. One can generally estimate F by measuring the recombina¬ 
tion emission induced by it [22, 23]. Although uncertainties connected with this quantity 
are large, their relevance is largely beyond the scope of this paper, since metagalactic flux 
is the dominant ionization source only at large heights above the plane of the disk and 
at large galactocentric radii, i.e. in the halo gas. We briefly mention that an interesting 
candidate for ionizing source present in all galaxies could be the massive decaying neu¬ 
trino, which not only solves the problem of the ionization, but also accounts for the dark 
mass present in every spiral galaxy. There exist a large number of alternative sources of 
intergalactic ionization: population HI stars, elliptical galaxies, dwarf galaxies, starbursts, 
quasars concealed by dust obscuration in intervening galaxies, reflection-dominated hard 
x-ray sources, accreting black holes and supernova-driven winds from early galaxies ([1], 
and references therein), or non-radiative sources like large-scale intergalactic shocks. Most 
of them, fortunately, are not operational in late epochs of galactic evolution, and their 
respective influences will be considered in the extension of the present work. 


The signihcance of these results lies mainly in exciting possibility of direct testing of 
models through high-sensitivity observations of recombination lines [24]. One of problems 
facing any such observations, especially in cases of small inclination angles, is strong signal 
originating in the fluorescing haloes interfering with light from the disk. We roughly 
compare two intensities using the estimate of local {z = 0) intensity of Ha halo emission 



where X is the value of the numerically solved integral J n^ds = 2.75T2' ^la cm ® pc [26] 
and expression for Ue is taken from the rhs of the eq. (1) using hrst two terms. A^hi is 
the neutral hydrogen column density along the line of sight through the halo only. In the 
Figure 1. we plot the dependence of /q, that is interstellar Ha intensity in rayleighs (1 
R=10®/47r photons cm“^ s“^ sr“^), on the Galactocentric distance, r. Although poorly 
known, A^hi is expected to be ~ 10^® cm“^, characteristic of the average Lya absorption 
originating in haloes of normal galaxies at low redshift [25]. Therefore, the signal from the 
disk should be strong enough to be analyzed in most of the typical cases. Of course, the 
kinematic properties of the emission lines originating in halo and disk will also be different, 
helping practical discrimination between the two. 

We note that Reynolds et al. [27] obtain for /q, at 1^1 ~ 1 kpc the following value: 
47r J ~ 2 X 10® photons cm ^ s where J is the incident Lyman continuum flux, that 
can be compared with our estimate depicted in the Figure 1. The project that is currently 
operating - Wisconsin Ha Mapper (WHAM) [28] will provide the data on the diffuse 
ionized hydrogen through the optical Ha line. In a subsequent work, we intend to include 
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all terms in the eq. (1), in order to get complete picture of expected Ha emission, by 
solving the integral f n'^ds. 
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Figure 1. The interstellar Ha intensity in rayleighs is plotted against the Galactocentric 
distance r in kiloparsecs. Different values for the height above the Galactic plane, z (where 
z goes from the disk to 1 kpc) are taken into account. 
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